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Recently a new propagating mode of coupled charge and spin oscillations was predicted in a two 
dimensional electron gas with a sufficiently strong Rashba interaction. We show that Coulomb inter- 
actions qualitatively modifies the spectrum and increases the characteristic wavelength of this mode 
by orders of magnitude, but does not suppress it. An absorption experiment that can conclusively 
detect the presence or absence of such a propagating mode is proposed. 
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In the field of spintronics, semiconductor devices uti- 
lizing the spin degree of freedom enjoy special attention 
due to prospects of their ready integration with exist- 
ing solid-state technologies "l^ . In addition, the compar- 
atively strong spin-orbit interaction in semiconducting 
materials allows for the possibility of spin manipulation 
by electric fields. This ameliorates the challenge of cre- 
ating local and variable magnetic fields in the nanoscale 
regime to otherwise operate arrays of spintronic devices. 
The study of the interplay of spin and charge motion in 
semiconductors (see review 0) dates back to research on 
optical spin orientation j^. Subsequently, spin-galvanic 
effects were discovered, in which uniform spin polariza- 
tion created by optical orientation or other techniques 
generates electric currents (see review 0). An inverse ef- 
fect of inducing non-equilibrium spin density by current 
was suggested theoretically §] and studied by several 
groups 0, H, 0, 1^- The spin-Hall effect, predicted in 
Ref. Illl is another phenomenon in which electric current 
creates a spin response by generating a pure spin current 
in an orthohonal direction. A considerable discussion 
(see review fl^ of the nature of such a spin current is un- 
derway following two recent puclications 0, Q| on the 
subject. 

In this paper we consider a specific case of a Rashba 
two-dimensional electron gas (2DEG) system in the dif- 
fusive regime. With regard to potential applications, this 
regime is experimentally much less challenging than the 
ballistic regime, required in many other spintronic pro- 
posals. Moreover, in the 2DEG the spin-galvanic and 
spin-Hall effects are clearly distinguishable by the spin 
direction involved: the former is relevant for the in-plane 
spin density, while the latter creates a current of out- 
of-plane spin. When an alternating current is passed 
through a 2DEG, additional features may appear. Re- 
cently, a new propagating mode of coupled spin and 
charge oscillations was predicted in 2DEG structures 
with a sufhciently strong Rashba interaction ^3- ^ 
low frequency mode, existing in the diffusive regime, is 
of major potential importance in semiconductor physics. 



However, the derivation in Ref. 113 treated electrons as 
neutral fermions. Such a treatment can be justified in 
the stationary cases 0, 0| , but not in the case of the 
time dependent electron 5n{r,t) and associated charge 
p — eSn density deviations in the propagating mode. In 
bulk conductors the presence of the Coulomb interaction 
prevents the existence of low frequency modes by impos- 
ing the quasi-neutrality condition, p = 0. An important 
exception are modes propagating in large external mag- 
netic fields e.g., helicons. Since screening in two 
dimensions is weaker than in three dimensions, the quasi- 
neutrality condition is relaxed. But the uncompensated 
charges will create an induced electric field Eind which, 
in turn, will modify the motion of the system. A proper 
self-consistent account of the induced electric fields in the 
proposed mode is the main goal of the present paper. 

Following Refs. EE E 113 we consider a 2DEG with 
the single-electron Hamiltonian 

n2 

H = h az\(j X p] 

2m 

Here m is the effective mass, a is the Rashba interaction 
constant, £ is a unit vector perpendicular to the plane of 
the 2DEG, and a is the vector of Pauli matrices. 

For temperatures much lower than the Fermi energy, 
assuming an isotropic momentum relaxation time Tp, the 
linearized diffusive transport equations for two dimen- 
sional charge and spin densities 6n and S, respectively 
have the following form 0, 0, E| : 

Sn = Ddk{dk5n-veEk) + 2TeikdiSk 
I An Q. r 

V Ts Ts 2 

= D{dkfS. - ^^f^dk~Sk - — (1) 
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Here i,k € {x,y}, v = m/{TTh^) is the density of states, 
and D = VpTp/2 is the diffusion coefficient. The param- 
eter r — 2a(eQTp)^/(l + 4:{eaTp)^) and the spin relax- 
ation time Ts = Tp(l + 4(eQ,Tp)^)/2(eQrp)^ are expressed 
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through the spin-orbit sphtting = app , where vp and 
pf are the Fermi velocity and wave vector, respectively 
|l7|. The renormahzation S = S — lygfig'B, where g is 
the Lande's g factor, accounts for the Zeeman action of 
the magnetic field. It is only significant in the presence 
of a large externally applied magnetic field and can be 
neglected here. 

The electric field E = ^ext + ^ind in Eq. is a sum 
of the externally applied field driving the system, and 
the field induced by charges and currents in the 2DEG. 
First, we explore the effect of the induced electric field 
on the spectrum of the propagating mode. The speed of 
this mode is much smaller than the speed of light [l5| . 
hence the induced field is quasistatic: Ei„d = —V4>ind + 
A^ind/c with potentials determined by the instantaneous 
values of electric charges and currents. The first term 
describes the Coulomb repulsion of induced charges and 
the second term accounts for the self-induction of the 
induced currents. The latter is small compared to the 
former, and will be neglected. As a result 



Er'^ir) = J T,{r - v')e5n{v')d?r' 



(2) 



where r is the in-plane coordinate and the form of the 
Coulomb kernel Ti depends on the geometry of the exper- 
iment and the presence or absence of the gate electrode, 
as will be discussed below. 

The propagating mode is a periodic solution of Eqs. 
with a density component ^ne*'^'^~'"* and a single nonzero 
spin component ^^^(fc, w)e*'^^~*"* (we choose i||k) il5| . 
The induced electric field of this mode is directed along 
k. Introducing the spin-diffusion length Is — \/Dts, and 
defining dimensionless quantities q = klg, qs = ve^ls, 
and 7 — Ty'Ts/D, we perform Fourier transforms on 
Eqs. with the result, 



iuj-4>i/Ts 2i-fq/Ts 
i-^4>il2qTs iui~(l)2lTs 
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where F is the Fourier transform of the kernel of Eq. Q 
F{k) = J Tx{'r)e~^'^^dxdy, and o-q = ye^D. The spec- 
trum of the mode is then given by 



Wl,2 = I — 

Ts 



-i{<t>i + 02) ± \/W¥l 



(4) 



Coulomb interaction enters the spectrum through the 
renormahzation factor R — l + iFqe/q with large param- 
eter qB ^ I, refiecting the strength of the Coulomb re- 
pulsion. For the case considered in Ref.fl^L D = lOcm^/s, 
Ts = 1 ns, Tp = 1 ps, m/me = 0.07 we estimate qs ~ 450. 

For small wave vectors both wi,2(fc) are imaginary, cor- 
responding to two purely dissipative modes of the sys- 
tem: the charge diffusion and the spin relaxation. As k 



is increased, wi_2(fc) acquire a real part when the expres- 
sion under the square root in Eq. becomes positive, 
and the mode can become weakly damped if the condi- 
tion Re{ijj) / Im{ijj) 1 is satisfied 0|. The existence of 
such a regime depends on the value of 7, and the proper- 
ties of the function F(k). First, consider a single 2DEG 
layer of thickness . Since the diffusive transport equa- 
tions are only valid for large wavelength 27r/fc 3> I, where 
I ^ vpTp^ dz IS the electron mean free path, we can ex- 
pand F{k) in powers of small kdz and obtain: 



Fs{k) 



~2m - 



Second, consider a 2DEG with a gate electrode located 
at a distance dg > d^. Here the kernel becomes smaller 
due to partial screening by image charges. Expansion in 
kdz ; kdg <C 1 gives 



-27ri I dgk 



d,k 



Sketches of the uji,2{k) functions are shown in Fig. ^ 
The spectra in the single-layer and gated cases have the 
same overall shape. Due to Coulomb repulsion, the real 




FIG. 1: Schematic plots of the real (dashed line) and imagi- 
nary (solid line) parts of aji,2(fc). A - with Coulomb repulsion, 
B - neutral carriers, C - ratio Re{ijj) / Im{Lj) corresponding to 
plot A, D - same ratio for plot B. 
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part of a;i_2(fc) exists in a finite interval [k-,k+], in con- 
trast to the case of neutral carriers where it exists for 
k > k":. The mode propagates with least dissipation for 
the wave vector /cmax at which |i?e(ix')//TO(tt')| reaches a 
maximum. One can prove that, in the physically relevant 
regime, qb 3> 7, the Coulomb renormalization in a single 
layer 2DEG, Rs = 1 + 2TrqB/q « 2iTqB/q ^ 1, is large in 
the relevant interval of wave vectors. In the gated 2DEG 
Rg = 1 + 2'nqB^dg — dz/i)/ls ^ 1 also holds for typical 
systems. In the limit of large Coulomb renormalization, 
our calculation gives the following values for the charac- 
teristic parameters of the curves shown in Fig. A) . 




Im{uj±) = 

-Re(Winax) ^ 



7^ + 1 ±7^72 + 1 



Re{Lo) 



(5) 



Im{uj) 



For neutral carriers the characteristic parameters (see 
Fig-IHB)) are 




FIG. 2: Numbered lines represent absorption curves for rep- 
resentative wave vectors fci < A;2 < . . . < fee. Thick solid lines 
correspond to = 0. As fc is increased, the position of the 
maximum of absorption moves to higher frequency until the 
maximum disappears. A - case of 7 > 1, B - case of 7 < 1. 
Inset of B - case of 7 = 0. 
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Eqs. (|3J) and © show that the Coulomb interaction in- 
creases the wavelength of the propagating mode. At the 
point of least damping the wavelengths are 27r/fc^a^x = 
An^qBls and 2n/k^^^ = 2^^ls^/R^ in the single and gated 
cases, respectively. Both wavelengths are much larger 
than for the neutral carrier case, 27r/k'^^^ ~ Is- The 
increased wavelength allows the system to reduce the 
Coulomb repulsion. Accordingly, a smaller increase is ob- 
served in the partially screened gated case (V^ < 2nqB)- 
Still, in both cases the wavelength can become compa- 
rable to the size of a typical 2DEG sample, at which 
point the boundary conditions, not considered here, will 
start to play an important role. The condition for the 
mode to be weakly damped, however, is unchanged by 
the Coulomb repulsion and remains 7^1 ^2111. Note 
that the diffusive regime condition, 27r/fc ^ I, is well 
satisfied in the large R case. 

A spatially periodic external electric field can excite 
the propagating mode and one expects that in the weakly 
damped regime a maximum of energy absorption at the 
frequency u = Re{uji^2{k)) will be a signature of the 
mode's presence. Experimentally, the spatial periodicity 



of the driving field can be created by an array of gates, as 
in experiments on plasmon resonance [20j . or, optically, 
using crossed light beams. The average dissipated power 
equals Q — E^xtR^ij) /2, where the Fourier component 
of the current j{k,uj) can be found from the continuity 
equation eSn + divj — as j — e{uj/k)6n. Finding Sn 
from Eq. JSJ one obtains: 



Qiw) 



w'^iw'^ + A) 



w 



A = 
C = 



7'(0i +202)) 
,2 
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(7) 



hi - 27^01 
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Representative absorption curves are shown in Fig. 12 For 
7 > 1 the peak is indeed observed for k values support- 
ing the propagating mode. However, Q{uj) also exhibits a 
maximum for small wave vectors where a;i^2(fc) are purely 
imaginary. The maximum is present even in the 7 <IC 1 
case. The exact condition for the existence of the maxi- 
mum is A > B or, equivalently, k < k^, where fc* is deter- 
mined by the equation 72(30i(fc*) -|- 4>2{k*)) = 0i(fc*)^. 
All maxima exist in the diffusive regime, 21: jk^ ^ I. 

To understand the presence of an absorption maxi- 
mum in a system with two purely diffusive modes, note 
that the spin-orbit interaction renormalizes the DC con- 
ductivity to a = (1-1- 27^) (To- At finite frequencies the 
conductivity decreases and reaches its bare value tro at 
1/ts <C <C 1/Tp. Accordingly, the fc = curves in 
Fig.|21are bell-shaped with a maximum at w = 0. But at 
any fc > the DC current is required to vanish. As a re- 
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suit, for fc ^ the curve Q{oj) starts from zero and then 
rapidly approaches the Q{u})\k=o curve, thus creating a 
pronouncedly asymmetric maximum, which looks more 
like a dip at w < Wmax(fc) (see Fig.|2lB)). In the absence 
of spin-orbit interaction the fc = curve becomes a hor- 
izontal line (inset in Fig. |2lB)). Thus an observation of 
the bell-shaped Q(w)|fc=o curve can serve as a proof of 
sufficient experimental sensitivity, allowing detection of 
the spin-orbit interaction. 

The maxima observed at 7 > 1 are much more sym- 
metric and significantly extend above the Q{uj)\k=o curve 
as shown in Fig. [SJA). Notably, the shape of (5(w) con- 
tinuously evolves with fc, showing no singularities at the 
boundaries k± of the interval of the propagating mode's 
existence. The distinction between the systems with 
and without the propagating mode has to be made by 
matching the experimentally measured Q{uj,k) familiy 
of curves to Figs.E^A) orE^B). 

To conclude, the propagating mode can be viewed as a 
periodic involvement of direct and inverse spin-galvanic 
effects and thus called a "spin-galvanic" mode. In this 
mode the spin-orbit interaction plays the roles of both 
inertia and a restoring force. Oscillations occur between 
two states with the same spin-orbit energy: one with 
non-equilibrium spin density and charge current, and the 
other with non-equilibrium spatial charges and spin cur- 
rent. The frequency of the mode is of the order to < 
1/ts <C 1/Tp. The spin-galvanic mode has no connection 
with the two-dimensional plasmon modes [2^ Esf . Al- 
though the spectrum of these plasmons w ~ Vk does not 
have a gap, they become damped and unobservable for 
Lu < l/Tp. The inertial element of plasmon oscillations 
is the mechanical inertia of electrons, thus they cannot 
survive in the diffusive regime. 

We have shown that the Coulomb repulsion of elec- 
trons modifies the properties of the mode, but does not 
suppress it. The low-frequency spin-orbit mode still ex- 
ists for r ^ -y/ D/ts, similar to the way the low frequency 
modes exist in magnetic field However, the Coulomb 
interaction modifies mode's spectrum and significantly 
increases its wavelength both in single layer and gated 
2DEGs. This has important implications for the experi- 
mental observation of the mode. On the one hand, spa- 
tially periodic electric field with larger wavelength can be 
created by purely optical methods, without the need to 
deposit periodic arrays of gate electrodes. On the other 
hand, it has to be checked that the wavelength of the 
mode does not approach the size of the sample. Theo- 
retical description of the latter situation will require a 
derivation of appropriate boundary conditions and a re- 
calculation of the resonant frequency. 

The original paper suggested that a narrow charge 
packet can be used to probe all wave vectors in the in- 
terval of the mode's existence. The evolution of such a 
packet would reflect the presence or the absence of the 
mode. Since an experimental method of injecting the 



excess charge into a 2DEG is not obvious, and since the 
domain of the spin-galvanic mode's existence is narrowed 
to a finite interval of wave vectors by the Coulomb inter- 
action, it may be beneficial to probe the mode at a given 
wave vector, as suggested in this Letter. We have shown 
that such an experiment is capable of conclusively prov- 
ing either the presence or the absence of the mode by 
observing the frequency dependence of the absorption of 
high frequency electromagnetic field, spatially periodic in 
the plane of 2DEG. 
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